Joint anisotropy and source count constraints on the contribution of blazars to the 

diffuse gamma-ray background 
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We place new constraints on the contribution ol blazars to the large-scale isotropic gamma-ray 
background (IGRB) by jointly analyzing the measured source count distribution (log A'^-log S) of 
blazars and the measured intensity and anisotropy of the IGRB. We find that these measurements 
point to a consistent scenario in which unresolved blazars make < 30% of the IGRB intensity at 
1-10 GeV while accounting for the majority of the measured anisotropy in that energy band. These 
results indicate that the remaining fraction of the IGRB intensity is made by a component with a 
low level of intrinsic anisotropy. 
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Introduction. — The origin of the isotropic gamma-ray 
background (IGRB), the observed all-sky diffuse emis- 
sion at MeV to GcV energies, remains uncertain. Some 
or all of this emission is expected to arise from astro- 
physical sources (e.g., active galactic nuclei, blazars, star- 
forming galaxies, millisecond pulsars, galaxy clusters, 
cluster shocks, and cascades from ultra-high-energy cos- 
mic rays; see [l|) as well as possible exotic sources (e.g., 
dark matter annihilation or decay ) ■ The Fermi Large 
Area Telescope (LAT) Collaboration has provided a new 
measurement of the IGRB energy spectrum Q with im- 
proved accuracy, covering the broad energy range from 
200 MeV to ^ 100 GeV. However, the energy spectrum 
of the IGRB is reported to be consistent with a power 
law, and thus provides little clue to the origin of this 
emission in the form of spectral features. As a result, the 
contributions of individual source classes to the IGRB 
are poorly constrained, severely limiting our ability to 
search for signals of new physics or to place constraints 
on emission from exotic sources. 

One way to tackle the problem is through population 
studies of resolved sources. Gamma-ray source classes 
with resolved members, such as blazars and millisecond 
pulsars, are obvious candidate contributors to the IGRB 
via emission from their yet unresolved members. The 
source count distribution in flux (log A'^-log S, the number 
of sources, N, per unit flux, 5") of LAT-detected gamma- 
ray blazars has recently been studied [1] down to fluxes of 
S'loo ~ 10~^° cm~^s~^, where 6*100 denotes the individual 
source flux above 100 MeV. For the first time, the log A^- 
log S is found to be well-described by a broken power law, 
and the position of the break and the slope of the log A^- 
log S below and above the break have been measured. 

It is possible to estimate the contribution to the IGRB 
of blazars below the LAT point source sensitivity of 
^ 10^^" cm~^s~^ by extrapolating the measured logA^- 



logS* to zero flux. For the energy range 0.1-100 GeV, 
Ref. Q reports that unresolved point sources contribute 
22.5 ± 1.8% of the IGRB intensity measured by 0]. As 
blazars constitute the vast majority of LAT-detected 
sources, this is a good indicator of the expected unre- 
solved blazar contribution, as well as a firm upper limit 
for sources that follow the measured source count distri- 
bution. In the following we use the terms blazars and 
unresolved sources interchangeably. 

Another constraint, which has not yet been explored, 
is provided by the level of anisotropy of the IGRB. Re- 
cently the first measurement of the small-scale anisotropy 
of the IGRB has been made , while in the last few years 
predictions have been derived for the anisotropy of many 
gamma-ray source classes, including blazars and galaxy 
clusters Q, millisecond pulsars Q, star- forming galax- 
ies and dark matter annihilation and decay U These 
source classes often produce similar energy spectra but 
very different anisotropics, suggesting that anisotropy 
analysis could be a powerful tool for distinguishing pos- 
sible IGRB contributors. 

In this Letter, for the first time we use the observed 
anisotropy information to constrain the properties of the 
source classes contributing to the IGRB. We calculate 
the intensity and anisotropy produced by the unresolved 
members of a source population whose detected members 
follow a broken power-law log A^-log S (such as the LAT- 
detected blazars). We then compare the predictions for 
this class of models with the measured IGRB intensity Q 
and anisotropy to identify the range of logAf-logS* 
parameters that are consistent with both measurements. 

The source count distribution. — The log A^-log S of 
sources detected by the LAT is compatible with a broken 
power law 
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where A is the normaUzation, Sh is the fliix where the 
power law breaks, and a and /3 are the power-law slopes 
below and above the break, respectively. The fluxes S 
and S'b are implicitly normalized to 1 cm"^ s~^. 

The log A^- log S" of the Fermi LAT sources has been 
measured in several energy bands For consistency 
with the anisotropy measurements, which were per- 
formed only above 1 GeV, we use the best-fit parameters 
for the log A''- log iS in the energy band 1-10 GeV: A = 
(3.61±0.17) X 10-" cm2 s sr-^, S'b = (0.23±0.06) x 10-*^ 
cm-2 s-\ a = 1.52 ± 0.15 and P = 2.38 ± 0.15 H H^]. 
In the following, S denotes the integrated source flux in 
the 1-10 GeV energy band. 

From a given logA^-logS", the contribution to the 
IGRB intensity from the unresolved sources can be es- 
timated by integrating the distribution from the source 
detection threshold down to zero flux: 

where St is the flux sensitivity threshold for point source 
detection. Due to the energy-dependent angular resolu- 
tion of the LAT, the value of St for a given source depends 
on its spectral index. However, this spectral index bias 
is small in the 1-10 GeV range We thus derive an 
effective St for the value / = 5.5 x 10-* cm-^ si~^ re- 
ported in Table 6 of Ref. , which, for the best- fit log A^- 
log S parameters, corresponds to unresolved blazars con- 
tributing 14% of the measured IGRB flux in this energy 
range [111. We flnd St = 3.7 x 10-^° cm-^ s'^ for the 
1-10 GeV band. 

The angular power spectrum. — The Poisson term of the 
angular power spectrum of the sources, Cp, can be cal- 
culated from the log A^-log S. It takes the same value at 
all multipoles and is given by 

This formula gives Cp in the units appropriate for the an- 
gular power calculated from an intensity map, i.e., units 
of intensity^ x solid angle, where intensity is in units of 
the number of photons per area per time per solid angle. 
Evaluating Eq. |3] using the best-fit log A^-log S param- 
eters yields Cp,prcd = 8.0 x 10-^* (cm-^ s-^ sr-^j^sr. 
Under the assumption that the sources are point-like and 
unclustered on the angular scales of interest, the Poisson 
term is the only contribution to the angular power. 

The predicted power, Cp^pred, has an uncertainty due 
to the propagation of the uncertainties in the parame- 
ters of the log A^-log S function. We estimate this uncer- 
tainty from the logA^-logS in the range 0.1-100 GeV, 
where the parameters have smaller statistical uncertain- 
ties, and then rescale it to the 1-10 GeV range. The 
rescaling was done by computing the flux conversion fac- 
tor between the two different bands: k = Sx-^/Su-^ = 



{y~'^~^^ — x-^"*"^) / (v~^'^^ — u~^'^^), where 7 is the av- 
erage photon index of the sources, and x—y and u—v 
are the edges of the two energy bands. Then the Pois- 
son anisotropics in the two bands are simply related by 
C™ = K^Cp". From the full covariance matrix [l3| of 
logiV-logS parameters in the 0.1-100 GeV range, we 
obtain = 0.54 x lO-^^ (cm-^ s'^ si-^f sr, 

which, assuming 7 — 2.4, can be rescaled to 5Cp-^° = 
<^Cp,pred = 0.8 X 10"^* (cm"2 s"i sr-i)2 sr, so that 
Cp,pred = (8.0 ±0.8) X 10-18 (cm-2 s-i sr-i)2 sr. 

Ref. reports measurements of Cp, obtained by av- 
eraging the angular power spectrum coefficients C^ over 
the multipole range 155 < £ < 504, in the energy ranges 
1-2 GeV, 2-5 GeV, 5-10 GeV, and 10-50 GeV. Using 
the same analysis pipeline as Ref. Q , we have also calcu- 
lated the anisotropy for the 1-10 GeV energy band for the 
foreground- cleaned data, which yields Cp^data = (11.0 ± 
1.2) X 10^1* (cm^^ s^i sr^i)^ sr. This value can be di- 
rectly compared with the predicted value derived above. 
The two values are compatible at the 2a level so that un- 
resolved blazars can account for almost all of the observed 
anisotropy. We discuss this point further in the next sec- 
tion. The 2(7 upper limit on the non-blazar anisotropy 

is C|^u = (Cp.data - Cp.pred) + \J ^Cp Aiita. + ^^P ^picd = 

5.9 X IQ-^^ (cm-2 s-1 sr-i)^ sr. 

Using the best-fit log iV-log S, we also compare the pre- 
dicted Cp with the anisotropy measurements in the four 
energy bands used in Ref. (TableH.Inthis case we use 
the rescaling method described above to calculate both 
the predicted mean values and their uncertainties. The 
derived 2a upper limits on the level of residual anisotropy 
in each energy bin are reported in Table HI These limits 
can be used to constrain models of astrophysical or ex- 
otic source populations, based on their predicted level 
of anisotropy. We note that the uncertainties, and, ex- 
cept for the 1-10 GeV case, the central values for Cp^prod 
used to derive these limits rely on the rescaling method 
described above, and thus on the assumption of an aver- 
age index for the sources. However, we flnd that varying 
7 from 2.2 to 2.6 produces only a small change. 

Finally, as a technical remark, we emphasize that the 
use of the dimensionful intensity angular power, rather 

TABLE I. Measured angular power from the foreground- 
cleaned data in different energy bands, Cp,data [3i predicted 
power from unresolved blazars, Cp,pred (this work); and 2a 
upper limits on the residual anisotropy, Cpjj (this work). 
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FIG. 1. Constraints on blazar log N-\og S parameters (break 
flux, Sh, and faint-end slope, a) from tlie intensity and 
anisotropy of tiie IGRB. Regions in which blazars provide 
100% of the observed IGRB anisotropy and mean intensity 
in the 1-10 GeV energy band are shown; the widths of the 
regions indicate the 68% confldence intervals. Below these re- 
gions blazars overproduce the anisotropy and mean intensity. 
Labeled contours show the fraction of the blazar contribution 
to the IGRB intensity. The best-fit la parameter region from 
the Fermi source count analysis is marked, along with the 
best-fit 5'b H (dot-dashed line). 



than the dimensionless fluctuation angular power, eon- 
veniently avoids the need to treat contamination of the 
anisotropy measurement by possible residual Galactic 
diffuse emission or instrumental backgrounds. These 
backgrounds are, to good approximation, isotropic, or 
vary only on large angular scales, and thus their contri- 
bution to the intensity angular povifer spectrum appears 
only at multipoles far below the range used to measure 
the angular power reported in In the following, when 
discussing the IGRB intensity, /igrb, we use the mea- 
surement given in Q . 

Constraints on unresolved blazars. — We now explore 
more generally the parameter space of the log A^-log S 
function to determine the region that is compatible with 
the measured anisotropy, intensity, and source count 
data. We define the parameter space of the source count 
distribution by the position of the break flux, 5'b, and the 
faint-end slope, a, of the logA^-log^* function at fluxes 
below the break flux. We fix the normalization and slope 
of the log A^-log S at high fluxes, as the efficiency in de- 
tecting point sources at high fluxes is ~ 1, and thus these 
parameters are well-determined (i.e., potential biases in 
these parameters are small). For each point in the Sb- 
a parameter space we calculate the predicted /igrb and 



Cp from the corresponding log TV-log S function. 

In Fig. [T] we show the region of the log TV-log S pa- 
rameter space in which blazars contribute 100% of the 
IGRB intensity (light blue) and that in which they con- 
tribute 100% of the angular power (dark yellow) in the 
1-10 GeV energy band. The widths of these regions show 
the 68% (Icr) confldence level regions, reflecting the re- 
spective 1(7 uncertainties in the measured Cp and /igrb- 
Above the light blue region, blazars contribute less than 
100% of the measured IGRB intensity; below this region, 
blazars overproduce the IGRB intensity. Similarly, above 
the dark yellow region blazars do not contribute the en- 
tirety of the measured angular power, whereas below this 
region they overproduce the anisotropy. We emphasize 
that the constraint from the anisotropy measurement is 
stronger than that from the intensity measurement ex- 
cept for at very high values of a. 

As noted in the previous section, the predicted Cp^prod 
from the best-flt log A^- log 5 agrees with the measured 
Cp,data to within 2a. We now ask "how well do the pa- 
rameters of the log A'- log S function inferred from Cp^data 
agree with those found from the source count analysis?" 
The best-flt la region of Sh and a for the blazar logA^- 
logS* given in [1, |l3l overlaps well with the la region 
inferred from Cp^data- This is a non-trivial result, as the 
measured anisotropy and source count distribution are 
independent observables, determined from independent 
data analyses. Here the agreement is at the la level 
since the errors on the 1-10 GeV log A-log iS* parameters 
shown in the plot are taken directly from [J, [lOj and are 
larger than the rescaled ones used in the previous section. 

There is a region of parameter space in which blazars 
contribute 100% of the IGRB intensity without exceeding 
the measured Cp; however, this region has a high break 
flux {Sh ~ 10~*cm~^ s~^) which is incompatible with the 
break measured from the source count analysis. Such a 
high break flux can be robustly excluded, as it would lie 
in the flux range where the source detection efflciency 
is close to 1, and thus this kind of feature is unlikely 
to have been missed. Taking the measured value of the 
break flux as an upper limit, we find that the contribution 
from blazars in the region allowed by Cp^data cannot be 
more than ^30% of the IGRB mean intensity (see labeled 
contours in Fig. [1]), a value which is in agreement with 
the results of the source count analysis alone. 

To further demonstrate how anisotropy data can be 
a powerful tool for distinguishing between multiple sce- 
narios we test an alternative fit to the blazar log A'-log S 
obtained by Stecker & Venters Uj . A notable feature of 
this alternative fit is that it can account for ~60% of the 
IGRB mean intensity. We have calculated Cp from their 
logN-logS 0,[ll and, using a threshold of 3.7 x 10"^° 
cm~^ s~^ (the same used in the rest of our analysis), ob- 
tain Cp = (2.4 ± 0.4) X 10^1^ (cm-2 sr-i)^ sr (the 
error reported on this prediction being likely an overesti- 
mate since it neglects the covariancc of the parameters). 
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FIG. 2. Cumulative contribution of blazars in linear (top) 
and log (bottom) scale to the IGRB anisotropy (dashed) 
and intensity (solid) for the Fermi best-fit log A'^-log S (E > 
100 MeV) as a function of source flux. 

This value is a factor 2.4 larger than the measured 
value, in significant tension with the anisotropy data. 

Comparing the measured anisotropy of the IGRB and 
the predicted anisotropy from blazars leads to another 
important conclusion. Since, for the best-fit source count 
distribution, blazars already account for ~ 100% of the 
observed anisotropy and, in intensity units, Poisson an- 
gular power is additive, the remaining component (or 
components) making ~ 70% of the IGRB intensity must 
contribute a low level of anisotropy in order to not over- 
produce the observed angular power. Interestingly, this 
can be achieved quite naturally since some proposed con- 
tributors to the IGRB, such as star- forming galaxies [1], 
are expected to contribute negligibly to the anisotropy. 
On the other hand, this result implies strong constraints 
on source populations with large intrinsic anisotropy. 

We note that the anisotropy and intensity contribu- 
tions from a source population have different depen- 
dences on the source count distribution, and conse- 
quently they represent complementary observables which 
are sensitive to different source flux ranges. This is 
demonstrated in Fig. [21 which shows the cumulative con- 
tribution to the intensity and anisotropy above 100 MeV 
as a function of source flux for the Fermi LAT best- 
fit logiV-logS* parameters. From the relative flatness of 
the cumulative flux distribution below the threshold flux, 
it can be inferred that the IGRB intensity contribution 
from unresolved blazars has only a weak dependence on 
the effective flux sensitivity. The cumulative anisotropy 
distribution, however, falls off more quickly below the 
threshold flux, so the anisotropy from unresolved sources 
is more strongly dependent on the sensitivity limit, and 
improved point source sensitivity is thus likely to have a 
more notable impact on the measured IGRB anisotropy. 

Additional energy bands. — We briefly consider this 
analysis in other energy bands. The range above 10 GeV 



is currently not suitable since the error on the measured 
Cp is large and the logA^-logS' is not well-constrained. 
A natural extension is thus to include the low-energy 
range down to 100 MeV. However, spectral index bias 
is non- negligible in this energy band, and the predicted 
anisotropy is quite sensitive to the shape of the source 
detection efficiency function, which characterizes this 
bias. In addition, there is not yet a measurement of the 
anisotropy below 1 GeV. For these reasons, an analysis 
of lower and higher energy bands is left for future work. 

Conclusions. — We performed a joint analysis of the 
source count distribution of blazars and the measured 
anisotropy of the IGRB in the energy range 1-10 GeV, 
and find that a consistent picture emerges in which un- 
resolved blazars account for only 30% of the IGRB 
intensity but ^ 100% of the angular power. The sources 
contributing the remaining ~ 70% of the IGRB intensity 
are thus constrained to provide only a small contribution 
to the anisotropy. Viable models of sources contribut- 
ing to the IGRB must satisfy the upper limits on their 
anisotropy that we reported in the last column of Ta- 
ble m These results demonstrate the power of anisotropy 
information for constraining the origin of the IGRB. 
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